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The mean king problem is a quantum mechanical retrodiction problem, in which Alice has to 
name the outcome of an ideal measurement on a d-dimensional quantum system, made in one of 
(d+l) orthonormal bases, unknown to Alice at the time of the measurement. Alice has to make this 
retrodiction on the basis of the classical outcomes of a suitable control measurement including an 
entangled copy. We show that the existence of a strategy for Alice is equivalent to the existence of an 
overall joint probability distribution for (d+ 1) random variables, whose marginal pair distributions 
are fixed as the transition probability matrices of the given bases. In particular, for d = 2 the 
problem is decided by John Bell's classic inequality for three dichotomic variables. For mutually 
unbiased bases in any dimension Alice has a strategy, but for randomly chosen bases the probability 
for that goes rapidly to zero with increasing d. 

PACS numbers: 03.67.-a,02.50.Cw 



I. INTRODUCTION 

The mean king problem was first introduced by Vaid- 
man et al. [IJ, and has since received a lot of attention 
as a basic quantum mechanical retrodiction problem: In 
the story, physicist Alice faces the mean king, who asks 
her to prepare a quantum system, on which his men will 
perform a von Neumann measurement in one of a speci- 
fied set of orthonormal bases. Alice is not present during 
this measurement, and knows neither the basis chosen 
nor the result obtained. She is then allowed a final check 
on the system (typically including entangled records of 
the initial preparation), leaving her with some classical 
measurement values only. She is then told which basis 
was used and is asked to correctly name the values found 
by the king's men. 

In [l[ the system was a qubit, and the three bases 
involved were the eigenbases of the three Pauli matri- 
ces. These are a special case of mutually unbiased bases 
("MUBs") , which means that after preparing a basis state 
of any of these bases, the probability distributions in all 
other bases will be uniform. It was subsequently shown 
that a maximal set of (d + 1) mutually unbiased bases 
exist for a d-dimensional Hilbert space, whenever d is a 
power of a prime, but to decide the existence of such a 
set for other dimensions (e.g., d = 6) has proved to be 
very hard. A status report on this problem, which has 
also many repercussions on other problems in quantum 
information is to be found in 0|. 

Note, however, that the basic statement of the problem 
makes no reference to the MUB property of the bases cho- 
sen by the king. So, supposing the game is to be played 
in d = 6, why should the king make things difficult (if 
not impossible) for himself by trying to find first a set 
of mutually unbiased bases in that dimension? Why not 
just pick any bases at random, say? Moreover, as we 
will see below, in the mutually unbiased case, Alice has 
a very simple way to compute a safe strategy. Again, 
a really mean king might make things more difficult for 



Alice here. These remarks go against another intuition, 
which would seem to make the choice of unbiased bases 
the meanest option for the king: Indeed, if Alice were 
trying to just make a measurement in one of the king's 
bases, hoping to pick the right one, in all but the lucky 
case her results will be totally useless. If the game were 
to be played many times, and Alice's aim was to be right 
as often as possible, she could improve her guesses using 
the correlations between bases. This gain is nullified in 
the MUB case. But the problem is not set like that: We 
demand of a solution that Alice is right in every single 
run, and this is not made easier in the least by the exis- 
tence of some statistical correlations between the bases. 
In other words, the intuition that unbiased bases are an 
especially mean choice by the king is fallacious. 

We therefore drop the assumption of unbiasedness and 
ask, for any choice of finitely many bases by the king: 
Can Alice find a strategy, consisting of an initial entan- 
gled preparation and a suitable measurement on the joint 
system after the kings men are through with their part, 
such that she gets the right value with probability one? 

Not very much has been done about this problem with- 
out assuming mutual unbiasednes. Some special cases 
have been discussed in @, ; S @- However, the avail- 
able studies apparently remain incomplete even in the 
qubit case (d = 2). 



II. SUMMARY OF RESULTS 

To state our main results, let us fix some notation for 
the rest of the paper. The system Hilbert space on which 
the king's men make their measurement will be denoted 
by H, and has dimension d. The number of bases chosen 
will be k, and the bases themselves will be denoted by 
4>b(i), for b = 1, . . . , k and i — 1, . . . , d. An important 
property of a choice of bases is the space 1Z of hermi- 
tian operators spanned by all the \$>b{i))($b(i)\- This 
space describes how many density operators we can dis- 
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tinguish with measurements in the given bases. Since 
J2 i = 1 for any basis, we expect only (d— 1)~ 

dimensions giving new information for each basis, so to- 
gether with the identity we expect 1Z to be {k{d— 1) + 1)- 
dimensional. If this number is achieved, we will call the 
chosen basis set non- degenerate. Of course, dim 1Z can- 
not exceed the dimension d 2 of the space of all hermitian 
operators on H, so for k > (d+ 1) every choice of bases is 
degenerate in this sense. The interesting property in this 
case is that dim TZ = d 2 , i.e., that the set of bases is to- 
mographically complete. Of course, for k = (d+1), which 
is the standard case, non-degeneracy and tomographic 
completeness are the same property. 
For any pair of bases, the values 

Pbc(i,j) = \ m^mm 2 a) 

are the joint probabilities of a pair of d-valued random 
variables, each of which is uniformly distributed. We 
say that a collection of k bases admits a classical model, 
if these probabilities are marginals of some joint distri- 
bution of all k variables (each taking d values). Since 
this property only involves the absolute values of scalar 
products, not their phases, and therefore captures only 
a small part of the information about the relative posi- 
tion of the bases, it is perhaps rather unexpected that 
the existence of a classical model is very closely linked to 
the existence of Alice's strategy. This is described in the 
following Theorem, our main result: 

Theorem 1 Let {&b(i)} be a collection ofk orthonormal 
bases in a d dimensional Hilbert space. Then 

(1) if the bases are non-degenerate (in particular, k < 
(d+1)) and the bases admit a classical model, then Alice 
can find a safe strategy in the mean king's problem with 
these bases. 

(2) Conversely, if the set of bases is tomographically com- 
plete (in particular, k > (d+1)), and if Alice has a strat- 
egy, then the bases allow a classical model. 

(3) In the case (1), Alice's strategy may begin with a max- 
imally entangled state, and in the case (2) a pure initial 
state is necessarily maximally entangled. 

Before going into the proof, let us see what this Theo- 
rem says about some basic examples. 

A. Mutually unbiased bases 

By definition a set of k bases in d dimensions is mutu- 
ally unbiased if, with the notation from ([1]), we have 

Pbc(i,j) = SbcSij^ + (1 ~ <^)^2- (2) 

From this a classical model is obvious, namely k statis- 
tically independent uniformly distributed random vari- 
ables. In order to compute the dimension of the span 
of the |<&b(i))($b(i)|, let us take the (kd) x (fc<i)-matrix 




FIG. 1: Possible range of triples (p a b(l, 1),P(, C (1, l),p ca (l, 1)) . 
The range of triples admitting a classical model is described 
as tetrahedron inside this body, with the same corners. 



of Hilbert Schmidt scalar products (denned for operators 
A,B by (A\B) HS := tv(A*E)) of these vectors, which is 
just the expression J5]), interpreted as a matrix Mu,cj. Its 
rank is the dimension we are looking for, and easily com- 
puted as k(d — 1) + 1, by determining all eigenvalues of 
M. Hence MUBs are non-degenerate for all k < (d + 1), 
and for any number of MUBs the mean king can come 
up with, Alice has a strategy. This result was previously 
obtained by another method in 0|. 



B. Qubits 

Another interesting special case, discussed in 0], is 
d = 2, k = 3. Choosing a basis in d = 2 is the same as 
choosing a pair of antipodal points on the Bloch sphere. 
Three bases are tomographically complete iff these points 
do not lie in a plane. The existence of a classical model in 
this case is one of the ancestral problems of quantum in- 
formation theory, namely precisely the existence of such 
models for three dichotomic variables characterized by 
Bell's original three- variable inequality [7|. The joint dis- 
tribution (HJ belonging to two bases 6, c is characterized 
(for d = 2) by the single number pb c (l, 1). Fig. Q] shows 
the possible range of triples (p a6 ( 1 , 1 ) , p fcc ( 1 , 1 ) , p ca ( 1 , 1 ) ) . 
The range of triples admitting a classical model, and 
hence a safe strategy for Alice, is described by Bell's 
inequalities as the tetrahedron inside this body. If the 
bases are chosen independently and with unitarily invari- 
ant distribution (Haar measure), the probability for this 
subset is exactly 1/3. This can be computed analyti- 
cally by reducing it to a problem of three independent 
uniformly distributed vectors on the Bloch sphere. 
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III. PROOF OF MAIN RESULT 

In the first round, Alice chooses a Hilbert space fC and 
prepares a density operator p on Ti ®1C, The first system 
is left to the king's men, who perform their von Neu- 
mann measurement in one of the bases $h, leaving a 
state (|$f,(i))($ h (i)|(g)l)p(|$ h (i))($f ) (i)|(g)l), conditional 
on their measured result being i. Finally, Alice will make 
a measurement on Ti ® /C, with some outcomes x £ X. 
This is described by positive operators F x onW®K, with 
J2 X Fx — The precise nature of the outcomes is irrel- 
evant. All that counts is that the value x provides Alice 
with a rule what to answer, if the king discloses that basis 
b was used by his men. We can express this by introduc- 
ing a "guessing function", but we might just as well take 
the rule itself as the outcome (possibly grouping together 
some outcomes leading to the same guesses). Hence we 
choose the outcome set 

X = {l,...,d} k ={x:{l,...,k}^{l,...,d}} 

with the interpretation that x(b) is the answer Alice will 
give, if her measurement gave the value "x", and the 
King discloses b. The requirement that she is right every 
time is the basic equation for p and F x we have to solve: 

tr(p(|# 6 (i)X*&(*)l ® l)F x (\$ b (i))($ b (i)\ ® 1)) 

(3) 

i,x 

At this point we can make the first simplifications. 
Suppose, for example, that Alice has found a solution 
using a mixed state p, and that \& is some unit vector in 
the support of p, so that | \I/)(\I/ 1 < Xp for some positive 
A. Then after replacing p by |\1/)( V E'|, the zeros in ([3]) will 
still all be in the right places, and since we chose $ as a 
unit vector, we have found a solution with a pure initial 
state p = |*)(*|. Next, we write * as 

d 

* = (1 <g> S)n with n = \ aa ) (4) 

a=l 

fl is the maximally entangled vector, and S is an opera- 
tor, whose matrix elements in a suitable basis of fC are the 
vector components of 'J, normalized so that tx{S*S) = 1. 
Then in J3]) we can commute S past the projections act- 
ing only on the first factor, and simplify the expression 
by introducing the vectors 

$ 6 (i) - (|$ 6 (i)X*6(»)| <8 

= $ 6 (*) ® <$> b (i) en®n, (5) 

where the bar indicates componentwise complex conju- 
gation in the basis \a), in which fl takes the form J4|. 
Then the basic equation |(3]) becomes 

® S)*F X (1 <g> S)|S 6 (t)) = Xi.flr^fi) • (6) 



Now suppose 77 is a vector in the support of the operator 
in this bracket. Then substituting \r))(r)\ for the operator 
will still give zero, whenever i ^ x(b), and hence 

($&(t)|»7)=0 if i^x(b). (7) 

But also the scalar products for i = x{b) are essentially 
fixed: We have J2i \®b(i))($b{i)\ = I for any basis b, 
which translates to J2i $b(i) — ^ Yla ©• Therefore, we 
can sum ([7| over i, obtaining 

(h(i)\v) = (0|»7> *<,*(6). (8) 

We will call such vectors safe vectors for Alice (and the 
particular outcome x). 

A. Structure of safe vectors 

How many safe vectors can Alice find? A key role 
for answering this question is played by the space 1Z 
introduced in Sect. HU or, equivalently by its image 
1Z C Ti. ® Ti. under the identification of operators on 
Ti and elements oiTi®Ti: 

K =lin R {$ b (i)}, (9) 

its complex linear span 72. c, an d its orthogonal comple- 

■ — ± 

ment 72 . For every x £ X we arrive at the following 
alternative: It may happen that there is no vector 77 sat- 
isfying $8$ with (fi|»7) 7^ 0. Then all solutions of that 

equation are in 72 , which also means that such values 
x can never occur as a result of Alice's measurement. Al- 
ice's strategy will have to rely on the other cases, i.e., the 
subset of those x £ X, for which a non-trivial solution r\ 
of ([H| exists. To get a standard solution, we multiply r\ 
with a scalar so that (fl\rj) = 1. Moreover, we can ap- 
ply to T] the orthogonal projection to 72c, thus obtaining 
a solution which is uniquely determined, since all scalar 
products with vectors from this space are fixed. We note 
that since all its scalar products with the $b(i) are real, 
we can even conclude that r\ x £ 72. Hence whenever a 
non-zero solution exists for some x £ X, we can pick a 
unique solution r\ x , determined by the conditions 

($b{i)\Vx) = <5*,x(6) with r) x £H. (10) 

The fact that r\ x lies in this real-linear subspace means 
that the corresponding operator on Ti is hermitian, or, 
expressed in the standard basis that 

(aP\r) x ) = (pa\Vx). (11) 

It is clear that if Alice can find any safe vectors at all, she 
has some success at a unambiguous retrodiction game, in 
which she is allowed to pass, but has to be absolutely sure 
of her guess otherwise. As in the problem of "unambigu- 
ous discrimination" [8fl her aim would be to minimize the 
probability for pass moves. In the mean king problem, 
however, her success probability is required to be unity, 
which is the same as saying that J2x ^x — 1, and a guess 
x £ X is produced in every run. 
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B. Necessary conditions 

The Theorem states necessary conditions for the exis- 
tence of a strategy only in the tomographically complete 

case. Then 1Z = {0}, and the only choice Alice has is 
to pick safe vectors, which are multiples of the rj x as in 
(flQ|) . This fixes the operators in (J6j) to be 

(l®S)*F x (l®S)=p(x)\rj x )(ri x \, (12) 

with p(x) > 0. The values for x not allowing a non-zero 
safe vector can be subsumed by setting p(x) = 0. The 
overall normalization condition J2 X F x = 1 then reads 

N = Y,P(*)\v*)(rh\ = (l®S*S). (13) 

X 

Taking matrix elements of this equation in the standard 
basis and using the hermiticity (jTTJ) , we find 

(a/3\N\a'/3') = (f3a\N\f3'a'). (14) 

By lfl3|) this amounts to Saa'nppi = 5pf3>n aa i , where n 
is the matrix of S*S. With a = a' = 1 we find that 
S*S is also a multiple of the identity matrix. From the 
normalization condition trS*S = 1 this multiple must 
be l/d. Since S* S is just the reduced density operator 
of the restricted state, we have thus shown item (3) of 
the Theorem: In the tomographically complete case, the 
initial state of Alice must be maximally entangled. 

The connection with classical models is seen by taking 
the matrix elements of equation (fl3|) with other vectors. 
To begin with, let us consider the matrix element with 
f2. Then, since N — (1/(1)1, and (r] x \fl) = 1 whenever 
p(x) ^ 0, we get (n|JV|fi) = (l/d)<fi|fi> = 1 = E x p( x )- 
Hence the p(x) must indeed be a probability distribution 
on X, which is the same as the collection of all mea- 
surement outcomes. Furthermore, for any bases 6, c, and 
associated outcomes 

($ b (i)\N\$ c (j)) =$>(aO*i,*(6)* J - 1 x(c)- (15) 

X 

Clearly, the right hand side is exactly the marginal 
Pbc(i,j) of the probability distribution p with respect 
to the d-valued variables b and c. On the other hand, 
since N — (l/d) J, the left hand side evaluates to 

(l/d)($ b (i)\$ c (j)) = (l/d)\(Mmc(j))\ 2 , so p is ex- 
actly a classical model in the sense described in Sect. Hfl 
This completes the proof of the Theorem, part (2), and 
the corresponding statement in part (3). 

C. Sufficient conditions 

Let us now suppose, as in part (1) of the Theorem, 
that we are given k < (d + 1) bases. Then, for each x, 
(fit))) is an inhomogeneous linear system of equations for 
the vector rj x . Taking only the first d — 1 equations for 



each basis, plus one normalization equation (ft\i] x ) = 1 
eliminates the trivial dependencies between these equa- 
tions, so we have k(d — 1) + 1 equations for a vector in 
1Z . The condition of non-degeneracy described in Sect.lTTl 
is equivalent to saying that all these equations are non- 
singular, hence under the hypothesis of part (1) of the 
Theorem, r\ x exists for all x. 

Now suppose that a classical model exists in the form 
of a set of p(x) > such that the marginals lfT5|) are 
consistent with N = (l/ei)J, Note, however, that 1Z may 
now be a proper subspace, and the matrix elements with 
all $b(i) do not determine the operator TV completely. 
Nevertheless, we can set 

F x =dp{x)\r lx )(r lx \+F x (16) 

with F x > summing to the projection onto TZ . Then 
it is immediate that with a maximally entangled initial 
state, i.e., with the choice S = (l/y/d)l, the basic equa- 
tion ((6J is satisfied. This completes the proof of part (1) 
of the Theorem. 



IV. FINDING A STRATEGY NUMERICALLY 

Given the marginals of equation ([I]) , the existence of a 
classical model is a linear feasibility program in the p(x) . 
It can also be cast as a semidefinite program, namely 
to maximize J2 X p( x ) subject to the constraints p(x) > 
and ^2 x p{x)\rj x ){ri x \ < t/d. If the maximum turns 
out to be J2 x p( x ) = 1) we nave found the desired joint 
distribution. Otherwise, this is the probability for Alice 
to find an answer in the unambiguous retrodiction game 
described at the end of III. A. The following table lists the 
numerical results for low dimensions with independent 
Haar distributed bases, where ps is the probability that 
a safe strategy exists, E$ is the expected overall success 
probability for unambiguous retrodiction, and N is the 
sample size we used. 



d 


Ps 


E s 


logio N 


2 


.3334 


.6666 


7 


3 


.0013 


.398 


6 


4 





.34 
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Higher dimensions, with d d+1 variables and con- 
straints, are a serious challenge for PC based compu- 
tation. For d = 6, a strategy rarely exists, but one can 
first "debias" the bases with a gradient search minimizing 
J2ij a b Pab(hj) 2 - Instead of a semidefinite program one 
can then use the so-called EM-algorithm 0, [T0| to find a 
joint distribution, and this is typically successful for the 
debiased case, although convergence is rather slow. 
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